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We revisit previously developed analytic models for defect evolution and adapt them appropriately 
for the study of semilocal string networks. We thus confirm the expectation (based on numerical sim- 
ulations) that linear scaling evolution is the attractor solution for a broad range of model parameters. 
We discuss in detail the evolution of individual semilocal segments, focusing on the phenomenology 
of segment growth, and also provide a preliminary comparison with existing numerical simulations. 
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I. INTRODUCTION 

Cosmic strings are line-like concentrations of energy 
that can arise as topological defects in cosmological phase 
transitions in the early universe [l|, |21- Their notable 
popularity in the 80's was due to the belief that they 
could be responsible for structure formation, an alterna- 
tive scenario to inflationary perturbations l^-Q • At the 
end of the 90's, however, it was realised [y,0| that strings 
cannot be solely responsible for the observed structure, 
which was further confirmed by the direct observations 
of acoustic peaks in the CMB, in accordance with the 
inflationary prediction. This led to a shift of attention 
away from cosmic strings. 

However, recent developments in fundamental theory 
have led to a resurgence of interest in the subject [3, Q ■ 
First, it was realised that the production of cosmic strings 
is much more generic than previously thought in a wide 
class of cosmological models based on supergravity [lO] ■ 
In addition, new exciting developments in string phe- 
nomenology pTl - ll3l | allowed the study of cosmological in- 
flation from brane constructions, and this led to a generic 
picture of hybrid inflation ending in a phase transition 
that produces cosmic (super)strings [lj-[lg. 

These objects are much more complex than ordinary 
field theory strings and their properties depend sensi- 
tively on fundamental high-energy physics parameters of 
the underlying theory p7Hl9l |. Thus, not only it was 
realised that string-like defects are much more generic 
objects that could play a subdominant role in structure 
formation, but, also, their potential observability could 
open a window into physics at the highest energy scales. 

Over the past few years, there has been significant 
progress in understanding the properties, cosmological 
evolution and possible observational signatures of cos- 
mic superstrings (for a recent review see |20|). However, 
there remain certain types of well-motivated string ob- 
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jects for which their physics, and most notably their cos- 
mological evolution, remains relatively unexplored. For 
example, there is yet no detailed self-consistent model 
describing the evolution of cosmic strings carrying cur- 
rents, even though such strings have been predicted in 
situations where more than one cosmological phase tran- 
sition is considered [2l|, and they also arise naturally 
in SUSY models that predict strings [2z, 26\ . String cur- 
rents also appear in the effective 4D description of higher 
dimensional strings, as the string position in the inter- 
nal dimensions is described by worldsheet scalar fields 
giving rise to currents. Another example is semilocal 
strings [24j, arising in SUSY GUT inflationary models 
[25j and the corresponding D3/D7 brane inflation mod- 
els [2g. These models are a natural extension of usual 
inflationary models, in which the only extra ingredient 
needed is doubling a hypermultiplet. This simple change 
renders the properties of the strings very different; as 
we will discuss in the following section, these semilocal 
strings are non-topological, they have different stability 
properties [23: ^M ^^^ e^Q^ their CMB constraints are 
less stringent than those for usual cosmic strings |29| . 

More recently, monopole networks and hybrid net- 
works (where finite string segments have monopoles at- 
tached to their ends) have also been studied [30, W^ . 
Both in the gauge and in the global case, the hybrid 
networks tend to annihilate very quickly, usually in less 
than a Hubble time after formation, although they may 
survive for a longer period in particular circumstances. 
There is an outstanding hybrid-like case, where local 
strings attached to global monopoles. This, at least 
naively, is the starting point for the much more inter- 
esting semilocal case which we discuss here. 

String evolution is a difficult problem involving physics 
from different energy scales, which is often not fully un- 
derstood. There have been many different approaches 
for obtaining the properties of string networks, trying 
to negotiate between computational power and accu- 
racy, between numerical and analytic approaches. Nu- 
merical simulations approaches include field theoretic 
simulations, Nambu-Goto approximations^^ and the so- 
called unconnected segment toy models [a]- AH these 
approaches have their virtues and their limitations, and 



progress has been possible through combination of results 
from different fronts. Another complementary approach 
is that of analytic velocity-dependent one-scale (VOS) 
models [33,|33- These abandon a 'statistical physics' de- 
scription of the network in favour of a 'theromodynam- 
ics' one, characterising the network by a small number 
of physically meaningful macroscopic quantities such as 
correlation lengths and root-mean-squared (RMS) veloc- 
ities, whose evolution equations can be calculated from 
the microscopic equations. This approach has the ob- 
vious advantages of simplicity, tractability and applica- 
bility beyond the dynamical ranges accessible to simu- 
lations, but it comes at a cost — the averaging process 
(through which one deduces the macroscopic evolution 
equations from the microscopic ones) requires the intro- 
duction of a (relatively small) number of phenomenolog- 
ical parameters, whose values can only be inferred by 
using numerical simulations for calibration. 

In what follows we discuss such an analytic model for 
the evolution of semilocal string networks. We study the 
behaviour of the network as a whole, starting from the 
premise that it can be treated as a network of local strings 
attached to global monopoles. We also describe the evo- 
lution of individual semilocal segments, discussing under 
what conditions these segments can grow — a process that 
has been clearly identified in numerical simulations. Fi- 
nally, we provide a preliminary comparison between the 
analytic model and existing numerical simulations, focus- 
ing on the overall network properties and the dependence 
on the coupling parameter. A more detailed comparison, 
discussing in detail the distribution of the semilocal seg- 
ments and its time evolution, is left for a follow-up paper. 



II. SEMILOCAL STRINGS 



of section El reads 



S 



Semilocal strings [27|, |28|, |3J] were introduced as a min- 
imal extension of the Abelian Higgs model with two com- 
plex scalar fields instead of one, which lead to C/(l) flux- 
tube solutions even when the vacuum manifold is simply 
connected. The strings of this extended model will share 
properties with ordinary local U{1) strings, but they are 
not purely topological and will have different properties. 
For example, since they are not topological, they need not 
be infinite and can have ends. These ends are effectively 
global monopoles, which have some exotic properties by 
themselves |39| . 

It is thus clear that this semilocal string model is a field 
theoretical version of the type of object that we want to 
study in this paper (strings attached to monopoles). In 
section |Vl we will use numerical simulations performed 
using this field theoretical model to verify the results ob- 
tained in this paper, so we provide the reader with a 
short introduction to the semilocal string model in this 
section. 

The relevant action for the simplest semilocal string 
model, the one we will use in the numerical simulations 
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where $ — (0, %p). It can be easily seen that setting one 

of the two scalar fields to zero, we recover the Abelian 
Higgs model. We can therefore build from the analytical 
models applied to usual cosmic strings to tackle this new 
problem. 

The stability of the strings is not trivial, and it will de- 
pend on the value of the parameter /3 = '^tcaUr / "^^^a-ugc- 
for /3 < 1 the string is stable, for /3 > 1 it is unstable, 
and for /? = 1 it is neutrally stable. As we will see in 
section |Vl only low values of /3 will be of interest for 
the comparison, because otherwise the string network is 
either unstable or disappears very fast [40] • 

After a cosmological phase transition in such a model, 
it is expected that segments of semilocal strings will form. 
The cosmological evolution of a semilocal segment net- 
work will be quite different from the evolution of ordinary 
Abelian-Higgs strings [3J, |3^. The fact that semilocal 
strings have a different cosmological evolution is inter- 
esting because CMB predictions can be different [231 ^'^d 
can aid on inflationary model building [25|. 

The network evolution will depend on the interplay be- 
tween string dynamics and monopole dynamics. When 
the string segment ends, it must end in a cloud of gra- 
dient energy. Those string ends will behave like global 
monopoles that provide an interaction between strings 
independent of distance. Therefore, depending on the 
interplay between string dynamics and monopole dynam- 
ics, the segments can contract and eventually disappear, 
or they can grow to join a nearby segment and form a 
very long string f40| . 

We thus see that, at least to a first approximation, 
we can envisage these networks as being made of local 
strings attached to global monopoles, and as such the pre- 
viously developed analytic modelling techniques [30, 123 
should be applicable. This being said, it is also clear 
that these networks possess non-trivial dynamic al prop - 
erties, beyond those of standard hybrid networks (30l - [32 |. 
Specifically, the evolution of the string network will de- 
pend both on the string tension and on the dynamics 
of the gradient energy: the latter may be thought of as 
providing a long-range interaction between the strings. 
Nevertheless, we can hope to factor this in by modifying 
and extending these models so as to include the relevant 
effects. 



III. MONOPOLE NETWORK EVOLUTION 

A simple way to treat semilocal strings is by anal- 
ogy with the standard case of monopoles connected by 
strings, the difference being that the strings are gauged 
but the monopoles are global. We emphasise that in 
the present section we will describe the network by ex- 



plicitly modelling the dynamics and interactions of the 
monopoles. (In the following section we will follow the 
complementary approach and look at the evolution of in- 
dividual string segments.) This is particularly justified 
in the semilocal case since (as has been shown in previ- 
ous work [40]) it is indeed the monopoles that control the 
evolution of the network. 

The standard approach to analytical modelling of de- 
fect networks is to start from the microscopic equations 
of motion - e.g. the Nambu-Goto equations of motion for 
string defects - and perform a statistical averaging pro- 
cedure, under the assumption that the defects are ran- 
domly distributed at large enough scales. This allows 
one to construct a macroscopic energy evolution equa- 
tion (which can be traded for an equation for the net- 
work's characteristic lengthscale) and an equation for the 
RMS network velocities, which together describe the net- 
work at large-scales in a 'thermodynamical' sense. De- 
fect interactions are then added to these equations in 
a phenomenological way. In the case of cosmic strings, 
this procedure leads to the so-called VOS model [32, [Sj] , 
which has been well-tested against simulations and is 
used for predicting CMB signals of string networks. 

In the case of monopoles, one can follow an analogous 
procedure. The idea is to write an evolution equation for 
the monopole density, neglecting interactions, and then 
re-express it in terms of a characteristic lengthscale, L 
(the average inter- monopole distance). The effects of 
monopole forces and friction are then included in this 
equation (as well as in the relevant velocity equation) 
by adding extra phenomenological terms. More details 
can be found in [30, [3l| . The evolution equation for the 
characteristic monopole lengthscale reads as follows 
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where c* is a free parameter (to be calibrated by simula- 
tions) quantifying energy loss, and where we have defined 
a damping lengthscale. Id that includes the effects of ex- 
pansion (Hubble parameter H) and of friction (length- 
scale If) due to particle scattering 



The constants km and kg parameterise the monopole 
and string forces, and ris, ijm are the relevant symmetry 
breaking scales. Since we are mostly interested in late- 
time scaling solutions we will (unless otherwise stated) 
neglect the effect of friction due to particle scattering, 
which is only relevant in the early stages of the network's 
evolution. 

For most of this and the following section we shall 
consider standard expansion rates, corresponding to the 
parameter range < A < 1, so that A = 1/2 in the 
radiation-dominated era and A = 2/3 in the matter- 
dominated era. (Towards the end of the current sec- 
tion we will briefly comment on the limiting cases A = 
and A = 1.) The reason for explicitly studying radia- 
tion and matter era solutions is that recent observational 
constraints 29] show that semilocal string networks can- 
not be the dominant component of the universe's energy 
budget, but can only contribute a small fraction to it. 

Immediately we see a difference relative to the other 
hybrid networks. In [3l[ it was shown that for local 
strings attached to local monopoles the force due to the 
strings was always dominant when compared to that due 
to the monopoles, while for global strings attached to 
global monopoles the string force was subdominant at 
string formation but became dominant later in the net- 
work's evolution. In the present case the ratio is 
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and therefore the string force should always be subdomi- 
nant. This is in agreement with theoretical expectations 
and numerical simulations. 

Interestingly, this combination of terms means that the 
only attractor solution in an expanding universe (with 
a (X t^) is linear scaling 
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as in the case of plain global monopoles, and indeed we 
can closely follow the discussion in [30] • In particular, 
there are two possible branches of the scaling solution. 
First, there is an ultrarelativistic one with 
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The evolution equation for the RMS velocity v of the 
monopoles is 
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where the second term in the square bracket describes 
the force due to the strings, and for an expansion rate 
of the generic form a{t) ex t^, the Hubble parameter and 
Hubble distance are respectively 
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which only exists for relatively slow expansion rates (A < 
3/4) but is in principle allowed both on the radiation 
and matter eras. Second, a normal solution exists for 
any expansion rate, and its scaling parameters obey 
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As in the case of standard global monopoles there is a 
constraint on the velocities 
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FIG. 1: The time evolution of 7 (top panel) and v (bottom 
panel), for three different choices of initial conditions, with 
the parameter choices A = 1/2, fcm=l and c* = 0.23. 



which is trivial for A < 3/4 (that is, fo — > 1 is allowed), 
but restrictive for faster expansion rates. On the other 
hand, velocities will generically be significant: having 
fo — i- requires A — > 1. 

In what follows we briefly explore the parameter space 
of these scaling solutions. In Fig. [T]we show an example 
of the convergence of 7 and v towards the scaling solution 
in the radiation era, which for this particular choice of 
parameters is 70 — 0.0118, Vo — 0.0768. 

In Fig. [2 we illustrate the dependence of 7 on c*, with 
other parameters kept fixed: the slow growth is note- 
worthy. Since strings typically have lower energies than 
the monopoles, the energy loss parameter c* mostly re- 
flects the annihilation of monopoles and antimonopoles, 
with some contribution from radiative losses; loop forma- 
tion losses are expected to be much smaller. Therefore 
increasing c* corresponds to increasing the annihilation 
rate, and the typical monopole separation correspond- 
ingly increases. 

Figure [3] generalises the above and displays the depen- 




FIG. 2: The dependence of the scaling parameter 7 on c*, for 
A = 2/3 and fcm = l. 



dence of the asymptotic scaling values of 7 and v on the 
two parameters c+ and km, in this case for the matter 
epoch. As expected 7 also grows slowly with km, and v 
grows with both parameters. 

Finally, Fig. U quantifies the effect of the expansion 
rate, by displaying the asymptotic values of 7 and v as 
a function of c^ and A, for a fixed km = 1- Here the 
behaviour is more interesting. One sees that 7 has rapid 
and non-trivial variations for slow expansion rates, and 
inspection of the v plot reveals that this corresponds to 
the network being in the ultrarelativistic branch, and 
eventually (at moderate expansion rates) switching to 
the normal one. 

We note that it is also of interest to consider the case 
of Minkowski space (corresponding to A = and H = 0) 
but with a friction lengthscale proportional to the corre- 
lation length (say, for simplicity, £f ^ L). This should be 
an adequate description of some of the numerical simu- 
lations of semilocal strings done so far [401 ■ I^ ^^^^ case, 
linear scaling is still the attractor solution, although the 
scaling parameters now obey different relations, namely 
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It might be possible to test these relations numerically, in 
particular by comparing simulations done with different 
values of the parameter (3. 

On the other hand, in the opposite limit of fast expan- 
sion rate (A > 1, or in other words inflation) the linear 
scaling solution of Eq. [7] no longer exists. In this case 
the network is conformally stretched and approximately 
frozen, and its scaling solution is 
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These conformal stretching solutions are ubiquitous in 
the defects literature. 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 





0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



FIG. 4: The value of the asymptotic scahng parameters 7 
(top panel) and v (bottom panel) as a function of c* and the 
expansion rate parameter A, for a fixed km — 1. 



FIG. 3: The value of the asymptotic scaling parameters 7 
(top panel) and v (bottom panel) as a function of c* and km , 
for the matter-dominated era, A = 2/3. 



The analysis in this section may be simplistic in sev- 
eral ways, although the result of linear scaling is encour- 
aging, since it is seen in numerical simulations. The 
key issue that still needs to be considered is the role 
and consequences of the gradient energy. It may be 
thought of as providing a long-range interaction between 
the strings, but it is not clear what its dynamical effect 
on the monopoles will be. 

A possible way for describing the effects of the gradient 
energy is to treat semilocal networks as composed of two 
types of strings (standard strings and gradient strings). 
In this description an important issue is what happens 
at the string ends (monopole junctions) — there will be 
some matching condition that would need to be encoded 
in the evolution equations. 



IV. SEMILOCAL SEGMENT EVOLUTION 

Even if we can characterise the network by focusing on 
the evolution of the monopoles, it is still relevant to un- 
derstand the evolution of the individual string segments, 
although this may be possible only in an averaged sense. 
The evolution equations for string segments can straight- 
forwardly be derived in the context of the VOS model 
[111, yielding 
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where Is is the length of the segment under consideration, 
Vs its RMS velocity, k a parameter describing string cur- 
vature, and Id is now the string damping length, 



1 



= 2H 



1 



(16) 

























(a.,v^.)=(0.1,0.9) : 


































(a.,v^.)=(0.5,0.5) 


10"' 


- 








N 


N 




(a.,v .)=(0.9,0.1) , 














■s. 


^ 








10- 












s 


N N ■-. 


10-= 














N N ■ 
















S. N 


10 





































0.9 


/ 




- 


0.8 


/ 
/ 




- 


0.7 


- / 




- 


0.6 


/ 




- 


0.5 


t 






- 


0.4 


(a.,v^.)=(0.1,0.9) 


■ 


0.3 






{a.,v^,)={0.5,0.5) 


- 


0.2 






{a.,v^,)=(0.9,0.1) 


- 









10 
Vti 



FIG. 5: The time evolution of a = lajt (top panel) and Vs 
(bottom panel) for string segments with three different initial 
conditions, for a matter-dominated era and A; = 1. 



Again, we will neglect the damping effect of friction from 
particle scattering, and consider only that due to the 
Hubble expansion. At this point we are also neglecting 
energy losses in the segments due to radiation or loop for- 
mation; we will come back to these issues shortly. Nat- 
urally, in these circumstances all segments will quickly 
shrink, become ultrarelativistic and disappear, as illus- 
trated in Fig. [5] the segment size relative to the horizon, 
a — Is/t, and the segment velocity respectively follow 
a — ?> and Vs ^ 1. 

However, theoretical expectations, confirmed by nu- 
merical simulations, indicate that the dynamics of semilo- 
cal string segments is more interesting than that. While 
some segments do disappear, there is also a significant 
probability that segments merge and form longer seg- 
ments. This is due to the long range interactions of 
the global monopoles at the end of the segments, as 
has been discussed in Sect. H. In certain circumstances 
this process could even culminate in the disappearance 
of the monopoles, leaving behind a network of standard 
Abrikosov-Nielsen- Olesen strings . 

In what follows we discuss two possible phenomenolog- 



ical ways of accounting for this behaviour, by modifying 
the evolution equation for the segment size. A possible 
alternative would be to modify the forces in the veloc- 
ity equation, keeping in mind the notion that at large 
distances the force should be string-like (inversely pro- 
portional to distance) but at short distances it may be 
independent of distance. This alternative is left for sub- 
sequent work. 



A. Scale-dependent Behaviour 

A simple generalisation of the above equations would 



be 
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where L is the characteristic scale of the monopoles (and 
Vm their characteristic velocity) that has been previously 
discussed, and cr is a free parameter controlling the im- 
portance of the newly introduced term. We can in prin- 
ciple allow for a force due to the strings in the veloc- 
ity equation, but this is expected to be subdominant, so 
we will neglect it in what follows; a further simplifying 
assumption (justified from our analysis in the previous 
section) is to take I'm ~ 1 in Eq. (fTT)) . 

The new term was added on the purely phenomeno- 
logical reasoning that, to a first approximation, small 
segments should shrink and large ones should grow and 
merge. This can be intuited as a competition between 
two characteristic timescales. Each segment will have 
an annihilation timescale, and each monopole will have a 
characteristic timescale in which to find its (anti)partner, 
annihilate, and therefore produce a longer segment. The 
second process is expected to become relatively more 
likely as the segment size increases. 

One can search for scaling solutions of these equations, 
of the form a = Is/t = const., v^ = const.. Here we will 
assume that Vm — 1 and also that the overall density of 
the network has itself reached scaling, so that L — -ft with 
7 — const, as previously discussed. Again one finds two 
possible branches. There is an ultrarelativistic branch 
with 
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and also a normal one with 
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FIG. 6: The time evolution of a = Isjt (top panel) and Va 
(bottom panel) for string segments with different initial con- 
ditions in a matter-dominated era, with the parameter choices 
a = 1.2, k = \. 



Defining T = A [ivl - l) + 1, Eq. ^ has two real solu- 
tions, provided 



a > 4r7 ; 



if this is the case, then 
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Note that the negative sign solutions correspond to a ^ 
7. If (T is too small to fulfil the condition of Eq. ([22]). 
then the asymptotic behaviour is as in the standard case 
(a = 0). 

The above behaviour is illustrated in Fig. El for 
a = 1.2, a matter dominated era and the parameter 
choices c* — 0.23 (for the monopole network) and k — \. 
In this case the asymptotic values of a are respectively 
a = 2.787 and a = 0.7156, while the corresponding ve- 
locities are Vs = 0.2691 and Vs — 1. For some initial 
conditions one also observes a transient scaling solution. 
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FIG. 7: The scaling behaviour of a (top panel) and Vs (bottom 
panel) as a function of a, for segments of different initial sizes 
and velocities; the plots are for the matter-dominated era, 
with c* = 0.23 and k — 1. Also shown with dotted lines are 
the asymptotic behaviours in the hmits of large and small a, 
which are further discussed in the text. 



Here, therefore, the subsequent evolution of a segment, 
in particular which of the two branches of the scaling 
solution it follows, will depend on its initial conditions. 

One can also determine how the asymptotic scaling 
solution depends on cr, and under which conditions the 
above bifurcation occurs. This is shown in Fig. [T] 
In order to understand this behaviour one should note 
that when cr — > (while still fulfilling the condition of 
Eq. ^) we have 



1 



1 + A' 



that is, a — > 0. In the opposite limit, when ct — > cx) 
a . 1- A 



a 



(24) 



(25) 



1 - A ' ' 2Ao- ' 

so a — >■ 00 and Ws ^^ 0. For the parameters we have been 





FIG. 8: The scaling behaviour of a (top panel) and Vs (bot- 
tom panel) as a function of a and the expansion rate A, for 
a segment with initial parameters ai = 0.3 and Vai = 0.7; 
further parameter choices are c* = 0.23 and fc = 1. 



using (fc = 1 and A = 2/3) the latter equation yields 



1 

la 



(26) 



It is noteworthy that the bifurcation into the two dif- 
ferent branches (according to the initial conditions) only 
happens for an intermediate range of values of the param- 
eter a. Below this range the attractor solution is always 
the ultrarelativistic branch, while above it the attractor 
is the normal branch. 

Finally, Fig. [5] depicts the effect of different expansion 
rates on these attractor solutions. The asymptotic value 
of a increases significantly with A, particularly for high 
values of a. Recall that these solutions exist provided 
(J > AT J and that F itself depends on A. This threshold is 
more obvious for fast expansion rates, since this enhances 
the difference between the two solutions. Similarly for 
the velocities, the smaller the expansion rate the smaller 
will be the difference for values close to the threshold. 

In conclusion, with this phenomenological term, ini- 
tially small segments typically reach a scaling solution 



with smaller sizes and ultrarelativistic velocities, while 
large segments typically grow and reach nonrelativistic 
velocites. As a simple way to verify how robust (and 
model-independent) these results are, we will now con- 
sider an alternative form for this term. 



B. Balance Equation 

Here we will consider the following alternative form for 
the evolution equations 



dls 

dt 



= HL-vf^ 



dvs 
'dt' 



d— k. 



L 



(27) 



(28) 



In this we are assuming that the network of string seg- 
ments has a Brownian distribution, something that can 
be tested in numerical simulations. The new term ac- 
counts for the probability that different segments inter- 
sect, which depends both on the length/number density 
and velocity of the segments. In this case there are two 
new free parameters, d and ki. 

Before studying the scaling behaviour of this system, 
let us briefly describe how equation 1^7^ is derived. Con- 
sider a monopole network with characteristic length L 
and think of the monopoles as being connected by string 
segments of average length L. Now consider a string seg- 
ment of length Is < L moving with velocity Vs in this 
network. Assuming L has a Brownian structure, the 
probability that the segment Is will meet a segment of 
average length L in time St is: 



SP 



VsSt Is 






(29) 



Soon after the segments meet, a monopole from one seg- 
ment interacts with an anti-monopole from the other seg- 
ment and annihilates, so the segment Is joins with L in- 
creasing its length to ^^ -f L ~ L. Considering now a 
collection of segments of average length Is, so there are 
L"^//^ segments in a shell of volume L"^, the rate of change 
of network energy density ps due to such joining processes 
is: 



P. 



s,join 



= -d- 



sis fJ.{ls+ L) L^ Vs^J. 



L2 



L3 



/? 



d^, (30) 



where the coefficient d quantifies the efficiency of this 
joining process. This term, when translated to a rate 
of change of length Is gives rise precisely to the term 
with coefficient d in equation (j27p . The constant term 
fci describes the shrinking of the segment Is due to a 
monopole acceleration /„ = ki/ls giving rise to 



Ps, shrink 



/ifcl 



(31) 



We can now search for scaling solutions, as in the previ- 
ous case. Again we find two branches, an ultrarelativistic 
one 



Ws = 1 , Ct = 



fcl 



f-(l + A)' 



subject to the condition that 

d>7(l + A), 
and a normal one 

fci 



A(l-2t;2) + d^-l 



k 
2Aa 



which requires 



d > ^ [1 + A {2vl - 1)] 



Obviously scaling disappears if 



ki > d- 



.vJs 



(32) 



(33) 



(34) 



(35) 



(36) 



(37) 



A matter era example of the evolution of the segments 
in this model is shown in Fig. [HI Here large enough 
segments reach scaling (for this particular choice of pa- 
rameters this is a = 8.461, Vs = 0.08864), while for small 
segments (for which the sum of the d and ki terms is 
negative) there is no scaling at all, and the segments will 
shrink and disappear. Note that one can easily choose 
parameters such that a » 1, in other words the char- 
acteristic size of the segments becomes much larger than 
the horizon and the segments become rare; this can be 
seen in Fig. [TO] 

Finally, Fig. [TT] depicts the behaviour of the scaling 
solution (when it exists) as a function of several model 
parameters. As expected, a large c* (implying a large 
monopole separation) or a slow expansion rate make the 
existence of a scaling attractor more difficult, and can 
eliminate it altogether. In those circumstances all seg- 
ments would eventually disappear. This would then lead 
us to expect more and/or longer segments in the matter 
era than in the radiation one. 



V. COMPARING WITH SIMULATIONS 

It is vital to compare the results obtained from the an- 
alytical approach with some numerical simulations. The 
(simplest) model for semilocal strings that can provide us 
with some numerical simulations for comparison is the 
one introduced in section |lll As explained there, only 
low values of /3 are considered, because higher values ren- 
der the strings unstable or make the network disappear 
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FIG. 9: The time evolution a (top panel) and Vs (bottom 
panel) for three different segments in the matter era; further 
parameter choices are d — 0.02, fci = 0.5, c* — 0.23 and fe = 1. 



very fast. The actual simulations used in this work for 
a preliminary comparison to our analytic models are the 
(existing) ones described in [40|. We refer the interested 
reader to that reference, and only proceed to sketch that 
approach in this section (focusing on the aspects that are 
directly relevant for us). 

In a 512'^ box with periodic boundary conditions we 
discretize the equations of motion coming from Eq ([ij. 
As initial condition, we set up random velocities for the 
scalar fields in the semilocal model. After a transient 
time, the system forms a network that approaches a scal- 
ing regime. We are not interested in capturing physical 
details during the phase transition; instead, we are inter- 
ested in the long term properties of the network. This 
is why we only monitor the simulation to detect semilo- 
cal string segments once the system has relaxed to the 
scaling regime. 

Unlike ordinary cosmic strings, semilocal strings are 
not topological and that makes them more complicated to 
track on a lattice simulation. It can be the case that even 
though semilocal strings carry magnetic field strength. 
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FIG. 10: The asymptotic values of q (top panel) and Vs (bot- 
tom panel) in the matter-dominated era, for fci = 1, c* = 0.23 
and k — 1. 



FIG. 11: The asymptotic value of q as a function of d and 
c* (top panel) and of d and A (bottom panel) . In both cases 
the plots assumed fci = 1 and fc = 1; the top panel is for the 
matter era, while the bottom one has c* = 0.23. 



they have no zeroes of the scalar field in their core, so 
tracking zeros of scalar fields is not the best approach 
to detect the striiigs. The strategy used here is the one 
following [23, |40, SH : if the magnetic field strength at 
a given point exceeds a certain fraction of the maximum 
magnetic field of the corresponding Abelian Higgs cosmic 
string core, the point is considered to be part of a semilo- 
cal string. This would give us the 'volume' of the string, 
which we turn into a length dividing by the width of the 
corresponding Abelian Higgs string. Segments that are 
too short (typically shorter than a couple of times the 
width of the strings) are discarded. 

Therefore, we take snapshots of the simulated vol- 
ume at several time steps well within the scaling regime, 
search for lattice points which have significant magnetic 
field strength, compute how those points are collected 
into string segments, and calculate their length. It is 
the number density and length density of those strings 
that we use to compare to the predictions of the analytic 
models put forward in this paper. In what follows we do 
not directly compare the velocities in the model and sim- 



ulations, since reliable numerical measurements of these 
velocities are highly non-trivial and require the develop- 
ment of additional numerical algorithms, which must be 
left for future work. For analogous issues in the more 
standard case of Abelian-Higgs string networks, see [421; 
for the case of domain wall simulations with the Press- 
Ryden-Spergel algorithm see [43l |. 

We used results from several existing sets of simula- 
tions (of the order of 10 simulations for each value of the 
coupling /3 and the damping term d) in order to obtain 
basic statistics about the properties of the networks (de- 
tails of the simulations can be found in Table |l|. The 
fact that these are 'legacy' simulations (not done with 
the present purpose in mind) explains the difference in 
the number of simulations in each set, and also the fact 
that most of them have a constant damping term rather 
than a decaying one (which would more accurately de- 
scribe an expanding universe). Nevertheless, these are 
sufficient for our present purposes. 

As mentioned before, we are not focusing on the first 
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Simulations 


Time 


C 


M 


Is 


7m 


box=512^ 


150 


18821.8 ±378.8 


462.4 ±15.2 


0.56 


0.44 


S = 0.02 


200 


15462.9 ±314.9 


307.3 ± 14.0 


0.47 


0.40 


13 = 0.04 


250 


13426.6 ± 357.6 


222.4 ±9.3 


0.40 


0.34 


N = 9 


300 


11880.0 ±423.2 


167.8 ±9.4 


0.35 


0.31 


box=512^ 


150 


19811.3 ±529.2 


530.5 ± 12.8 


0.55 


0.42 


S = 0.02 


200 


16214.8 ± 293.6 


346.5 ± 14.4 


0.45 


0.36 


P = 0.06 


250 


14078.5 ± 110.1 


241.5 ±7.4 


0.39 


0.33 


iV = 6 


300 


12521.4 ± 142.6 


183.5 ±10.6 


0.35 


0.30 


box=512-'' 


150 


13535.3 ±362.4 


192.4 ±9.6 


0.66 


0.60 


S = 0.05 


200 


9888.8 ± 360.6 


103.1 ±10.2 


0.58 


0.55 


13 = 0.04 


250 


7624.3 ± 348.5 


56.4 ±6.7 


0.53 


0.53 


iV = 16 


300 


6254.0 ± 346.2 


36.9 ±5.9 


0.48 


0.51 


box=5123 


150 


12224.5 ± 345.4 


218.8 ± 13.1 


0.70 


0.57 


S = 0.05 


200 


9000.5 ± 482.5 


111.7±9.3 


0.61 


0.53 


13 = 0.06 


250 


7111.1 ±558.4 


62.4 ±8.5 


0.55 


0.52 


N = 23 


300 


5953.2 ± 572.4 


43.5 ± 7.4 


0.50 


0.49 


box=5123 


150 


14421.9 ±382.7 


244.5 ± 19.3 


0.64 


0.55 


S = 6/t 


200 


10166.8 ±398.6 


137.7 ±27.1 


0.57 


0.50 


13 = 0.04 


250 


7499.6 ±425.1 


93.1 ± 27.3 


0.54 


0.45 


iV = 16 


300 


5825.6 ±455.1 


74.7 ± 24.8 


0.51 


0.41 


box=512^ 


150 


12980.1 ± 196.2 


263.2 ± 14.1 


0.68 


0.53 


5 = 6/t 


200 


9155.0 ± 224.0 


124.7 ±11.3 


0.61 


0.51 


13 = 0.06 


250 


6813.2 ±336.9 


65.8 ± 7.9 


0.56 


0.51 


iV = 12 


300 


5413.9 ±327.6 


44.0 ± 8.2 


0.52 


0.48 



TABLE I: The measured values (with statistical errors only) 
of the total string length £ and monopole number TV for the 
various series of simulations further described in the text, and 
the inferred (approximate values of the string and monopole 
scaling parameters, 7s and 7m. 



time steps of the simulations, and the code was tuned to 
get into the scaling regime as fast as possible. Once in the 
scaling regime, the network was evolved using different 
(constant or time-decaying) damping terms. For each se- 
ries we follow the procedure described above, which gives 
us, for every simulation, a good approximation to the to- 
tal string length C and monopole number TV in the box. 
We proceed to average those two numbers over the vari- 
ous realisations, for several late-time boxes. The number 
of monopoles is obtained by counting the number of seg- 
ments and multiplying by two, since each segment has 
one monopole at each end. These will slightly overcount 
monopoles, since it does not take into account segments 
that have closed to form a loop. There are other system- 
atic errors as well in this 'number counting' approach. 
Therefore one should realise that the errors quoted are 
only statistical errors, and provide only a lower bound of 
the uncertainties in these simulations. 

We then translate the obtained string lengths and 
number of monopoles into a VOS-type lengthscale, which 
we calculate in two independent ways, either using the 
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FIG. 12; The evolution of the string scaling parameter 7s = 
Ls/t (top panel) and the monopole scaling parameter 7m = 
Lm/t (bottom panel), for four different series of constant- 
damping simulations. In each panel the four lines correspond 
to the following cases: 5 = 0.02, [3 = 0.04 (solid), 5 = 0.02, 
P = 0.06 (dashed), 5 = 0.05, /3 = 0.04 (dash-dotted), 5 = 
0.05, /3 = 0.06 (dotted). Further details are provided in the 
main text. 



strings or the monopoles: 

Ls 1 fV 



t 



1/3 



(38) 



Given the way they were calculated, Ls should be 
thought of as the typical inter-string distance (or perhaps 
the typical segment size), while Lm is a characteristic 
intcr-monopole distance. These are therefore not corre- 
lation lengths in the same strict sense as the term is used, 
for example, in Goto-Nambu string simulations. In par- 
ticular, the fractal distribution of the semilocal networks 
(and more specifically the assumption of a Brownian net- 
work) is an issue that warrants further study. Given the 
possible systematic uncertainties associated with the cal- 
culations of C and TV, we decided to present 7s and 7^ 
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FIG. 13: The evolution of the string and nionopole scaling pa- 
rameters for two different series of simulations with a damping 
coefficient S = 6/f. The solid and dashed lines respectively 
depict 7s and 7m for simulations with /3 = 0.04; the dash- 
dotted and dotted ones show the same quantities for /? = 0.06 
simulations. Further details are provided in the main text. 



without error bars. Figures [T^ and [T3] summarise our 
results. 

We start by noting the remarkable similarity between 
7s and 7™. The agreement is even more striking if one 
observes that, due to the previously mentioned possibility 
of an overcounting bias, the 7™ we obtain this way should 
be somewhat smaller than the correct one. Although we 
can't quantify the magnitude of this effect, inspection of 
animations of the simulations indicates that the number 
of loops is relatively small, and thus so should be the 
overcounting. Given the underlying uncertainties the two 
results can therefore be said to be in full agreement. 

We also find excellent agreement in the results of the 
simulations which only differ in the value of /3; there- 
fore the two available values of f3 are too close to each 
other to be distinguishable by this method. On the other 
hand, the effect of the different damping terms is quite 
obvious: a larger damping reduces the string length and 
the number of monopoles, thus leading to a larger 7. 
While keeping in mind the aforementioned uncertainties, 
we note that our results for the constant damping cases 
are consistent with the asymptotic value of the energy 
density in the network being inversely proportional to 
the damping, in other words p oc S''^. This is to be 
expected, for example if one considers the energy losses 
of a given length of string moving subject to constant 
damping. One may also expect the string and monopole 
velocities to decrease as the damping increases, but in 
this case the exact behaviour is not easy to determine, 
and one may also expect a noticeable effect of the param- 
eter /3. 

Finally, the absence of information on the velocities 
prevents us from providing a direct calibration of the pa- 



rameters of the analytic model. Nevertheless, we can use 
the results of the time-decaying damping case to attempt 
a naive back-of-the-envelope exercise. Referring to Eqs. 
(Pl IlOp we note that we have in principle three unknowns 
(the model parameters c^, and k„i, plus the velocity vq). 
We then note that the damping term 6 = 6/i (where in 
an expanding universe context this t of the simulations is 
effectively the conformal time) corresponds to a pysical- 
time expansion rate a (x t^ with A = 3/4. If we now use 
7 ~ 0.5 in Eqs. (|MT(I)l. we find the approximate relations 



C-ki^n 



3.2 



and 



Vq ~ 0.12fc„ 



(39) 



(40) 



These are plausible values. In the case of Goto-Nambu 
strings both the energy loss parameter c and the cur- 
vature paramater k are of order (but slightly smaller 
than) unity [33| . Here one may expect a slightly larger c* 
since the energy density of semilocal strings is somewhat 
smaller than that of Goto-Nambu ones: in the context 
of the analytic model, this would primarily be controlled 
by the parameter c^. As for our estimated velocity, this 
is somewhat lower than the ones typically encountered 
in other field theory defect simulations [42, |4j] , but one 
should also keep in mind that these previous measure- 
ments have almost always been performed in the matter 
or radiation eras — in other words, they have slower ex- 
pansion rates, and therefore an effectively smaller damp- 
ing term. 

In summary, although the above comparison is quite 
simplistic, the results are at least encouraging. As has 
already been pointed out, a more detailed comparison 
and a proper calibration of the analytic models must wait 
for the numerical implementation of a reliable method to 
measure defect velocities in semilocal string simulations. 



VI. CONCLUSIONS 

We have extended previously developed analytic mod- 
els for defect evolution, applying them to the study of 
semilocal string networks. We have shown that a linear 
scaling evolution (analogous to the well-known one for 
cosmic strings) is the attractor solution for a broad range 
of model parameters. Our results therefore provide sup- 
porting evidence for those of existing numerical simula- 
tions, where confirmation of the presence of this solution 
is restricted to the limited dynamical range available. 

We have discussed in some detail the evolution of in- 
dividual semilocal string segments, focusing on the inter- 
esting phenomenology of segment growth and considering 
two different possible scenarios thereof. As expected (and 
as seen in simulations) we find that small segments tend 
to shrink and decay while large ones can grow and merge 
with others. Which (if any) of the two scenarios we have 
considered is the valid one is an issue that can only be 
resolved with future and more detailed simulations. 



13 



Finally we have used the results of existing simulations 
in order to obtain a very preliminary comparison with our 
analytic results. This endeavour is quite limited by the 
fact that the available simulations have not measured the 
defect velocities, which play a key role in VOS-type mod- 
els. We found some encouraging results, but it is quite 
clear that a more detailed comparison, using expanding- 
universe simulations that explore a wider range of the 
space of model parameters (notably the coupling param- 
eter /3), will be needed in order to provide a meaningful 
calibration of the model. We shall address these issues 
in a subsequent publication. 
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